Phase relaxation of Faraday surface waves 
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Surface waves on a liquid air interface excited by a vertical 
vibration of a fluid layer (Faraday waves) are employed to 
investigate the phase relaxation of ideally ordered patterns. 
By means of a combined frequency-amplitude modulation of 
the excitation signal a periodic expansion and dilatation of a 
square wave pattern is generated, the dynamics of which is 
well described by a Debye relaxator. By comparison with the 
results of a linear theory it is shown that this practice allows 
a precise measurement of the phase diffusion constant. 
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Our understanding of spatio-temporal pattern forma- 
tion in non-equilibrium fluid systems has greatly bene- 
fitted 0| from recent quantitative experiments in com- 
bination with the development of new theoretical con- 
cepts. One of them is the so-called amplitude equation 
approach |^] , which is based on the linear instability of a 
homogeneous state and leads naturally to a classification 
of patterns in terms of characteristic wave numbers and 
frequencies. A different but equally universal description, 
the phase dynamics Q , applies to situations where a pe- 
riodic spatial pattern experiences long wavelength phase 
modulations. This approach, originally introduced in the 
context of thermal convection, has proven to be useful to 
understand the stability and the relaxation of periodic 
patterns, wave number selection, and defect dynamics. 
In many paradigmatic pattern forming systems such as 
thermal convection in a fluid layer heated from below 
(Rayleigh Benard convection, RBC) or the formation of 
azimuthal vortices in the gap between two rotating cylin- 
ders (Taylor Couette flow, TCF) the dominating wave 
number is dictated by the geometry and thus inconve- 
nient to be changed in a given experimental setup (for 
instance by a mechanical ramp of the layer thickness W). 

Faraday waves are surface waves on the interface be- 
tween two immiscible fluids, excited by a vertical vibra- 
tion of the container. Beyond a sufficiently large exci- 
tation amplitude the plane interface undergoes an insta- 
bility (Faraday instability) and standing surface waves 
appear, oscillating with a frequency one half of the drive. 
This type of parametric wave instability is attractive as 
the wavelength of the pattern is dispersion rather than 
geometry controlled. Just by varying the drive frequency 
the wave number can be tuned in a wide range. In that 
sense the Faraday setup is well suited for the study of 
phase dynamics. 

Nevertheless recent research activity in this field was 



mainly dedicated to the exploration of the processes un- 
derlying the selection of patterns with a fixed wavelength. 
Faraday ^ was the first to provide a quantitative study 
of this system, revealing that a sinusoidal vibration may 
induce a periodic array of squares. Later on, more com- 
plicated patterns with up to a 12-fold rotational symme- 
try (quasi-periodic structures) have been observed Hf7|. 
Here the amplitude equation technique contributed con- 
siderably to unfold the governing spatio-temporal reso- 
nance mechanisms. Applied to a set of modes with 
different orientations but fixed wavelength, |kj| = k 
the resulting set of Landau equations lead to a semi- 
quantitative understanding Q of pattern selection in this 
system. Motivated by these advances the idea came up 
to apply more complicated drive signals composed of two 
or more commensurable frequencies . That way the si- 
multaneous excitation of distinct wavelengths gave rise 
to novel surface patterns in form of superlattices p-pT[ . 
Only recently the phase information carried by the par- 
ticipating modes was found to have a crucial influence on 
the visual appearance of the convection structures |p^| . 

In comparison to other classical pattern forming sys- 
tems such as RBC or TCF, the phase dynamics in the 
Faraday system is much less explored. In usual Faraday 
experiments the drive frequency (or frequency compo- 
sition including relative amplitudes) is held fixed while 
the overall drive amplitude is ramped in order to record 
the bifurcation sequence of appearing structures. To our 
knowledge none of the previous investigations used the 
excitation frequency uj as the primary control parame- 
ter rather than the drive amplitude a. That way it is 
particularly simple to impose phase perturbations on or- 
dered patterns and to study their relaxation dynamics. 
Moreover, doing phase dynamics on the Faraday system 
has the additional advantage of rather quick relaxation 
times, which in typical setups are one and two order of 
magnitude faster than for instance in RBC. 

The present paper reports a systematic investigation of 
phase relaxation on Faraday surface waves. Our study is 
focused on the relaxational dynamics of an ideal surface 
pattern with a square tesselation. By evaluating the re- 
laxation time of the pattern in response to small changes 
of the frequency, the phase diffusion coefficient has been 
measured. The experimental results are found to be in 
good agreement with the findings of the linear theory 
evaluated for a system of infinite lateral extension. 

The experimental setup consists of a black cylindrical 
container built out of anodizcd aluminium, and filled to a 
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height h of 4.2 mm with a silicone oil (kinematic viscosity 
v = 21.4 x 10~ 6 m 2 /s, density p = 949 kg/m 3 , surface 
tension a = 17.3 x 10~ 3 N/m). In order to study finite 
size effects, we used three different containers, with inner 
diameters L\ = 265 mm, L-2 — 185 mm, and = 125 
mm. A glass plate covering the container was used to pre- 
vent evaporation, pollution and temperature fluctuations 
of the liquid. Furthermore, to avoid uncontrolled changes 
of the viscosity, density, and surface tension of the liquid, 
all the measurements have been performed at a constant 
temperature of 30 ± 0.1°C. The Faraday waves were ex- 
cited by an electromagnetic shaker vibrating vertically 
with an acceleration allowing for simultaneous amplitude 
and frequency modulations in the form a(t) cos uj(t). The 
corresponding input signal was produced by a waveform 
generator via a D/A-converter. The instantaneous ac- 
celeration was measured by a piezoelectric sensor. In 
a preparatory experiment undertaken with a sinusoidal 
(i.e. unmodulated) drive a cos cot the critical acceleration 
amplitude cl c (lo) for the onset of the Faraday instabil- 
ity was determined by visual inspection of the interface 
while quasi-statically ramping a at fixed u> — 2nf (see 
Fig. [I]). Throughout the investigated frequency interval 
70 Hz< / < 110 Hz the surface patterns, which appear 
at a supercritical drive of less than about 1.1 x a c , always 
consisted of an ordered square wave pattern, which - af- 
ter some healing time - was free of defects (Fig. [j]b). In 
order to study the dynamics of phase-perturbed patterns 
we have carried out measurements of the average wave 
number k(t) of the Faraday pattern in response to small 
changes of the drive frequency uj(t) around a mean value 
luq. The w-modulation has been accomplished in two dif- 
ferent ways: (i) By discontinuous jumps (back and forth) 
between frequencies coq — Alo/2 and u>o+ Alo/2 with a rep- 
etition period T between 100-300s, sufficiently large for 
the pattern to relax, (ii) By a sinusoidal modulation of 
the drive frequency according to co(t) — loq + Alu sin (tit), 
with T = 1/F — 2TT/fl between 2 and 1000s. (within the 
frequency range of our study, the response time of the 
shaker to small changes of the drive frequency is less than 
10 ms and thus negligible). In both cases the vibration 
amplitude a(t) was co-modulated in such a way that the 
instantaneous supercritical drive e = a(t) / a c [u>(t)] — 1 re- 
mained constant (see Fig. |l|a). The bandwidth Alo of the 
modulation needed to be confined to a few Hz in order to 
avoid the occurrence of dislocation-type defects. Under 
those conditions the square pattern remained practically 
ideal without perturbations (see Fig. ^b), just expanding 
and contracting (breathing) in a homogeneous manner 
with the modulation period T. To obtain the temporal 
wave number dependence a full frame CCD camera sur- 
rounded by a set of 4 incandescent lamps was mounted 
some distance above the container. About 100 pictures 
of the light reflected from the surface were taken at con- 
secutive instances of maximum surface excursion from 
which the spatially averaged wave number k(t) was ex- 
tracted by evaluating the position of the principal peaks 
in a 2-dimensional FFT. 
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FIG. 1. (a) Symbols denote the critical acceleration ampli- 
tude a c as measured at the onset of the Faraday instability 
at constant drive frequency u = 2irf. The solid line indicates 
the theoretical prediction for the material parameters given. 
Along the dashed line the reduced drive amplitude is e = 9%. 
The phase relaxation experiments were performed by impos- 
ing simultaneous small changes of / and a, as indicated for an 
example by the solid bar. (b) and (c) show a square pattern 
and the associated 2-dimensional power spectrum as observed 
during the modulation experiments. 

The wave number k(t) followed the modulation in 
a relaxational manner. For the discontinuous modu- 
lation (i) this is directly apparent from Fig. ||. Here 
the relaxation time r has been derived by fitting the 
exponential decay of the data. In the type (ii) ex- 
periment k(t) oscillates around a mean value kg with 
an amplitude Ak rn and a temporal phase lag 6 (see 
Fig. [|). Introducing the complex wave number incre- 
ment Ak* = Re[Afc*] +iIm[Afc*] = Ak m e iS : its real 
and imaginary parts are plotted in Fig. |J as a func- 
tion of the modulation frequency F. The solid curves 
of this figure are fits of a linear Debye relaxator, where 
Afc* = Ak(Q = 0)/[l + iQt]. Fig. H shows the relaxation 
time r, as obtained from both types of experiments, as 
a function of the mean wave number fcoi the container 
diameter L and the reduced drive strength e. A depen- 
dence on the (small) modulation amplitude Aw could 
not be detected. The experimental data reveal a linear 
increase with the wave number ko and a proportionality 
to the square of the container size L. 
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FIG. 2. Temporal decay of the wave number k(t) in re- 
sponse to a step-like [type (i)] change of the drive frequency 
/. Also the drive amplitude a was changed to keep the re- 
duced amplitude e = a(t) /a c [ui(t)] — 1 constant. The decay 
of the wave number k(t) could be fitted by an exponential. 
(Parameters: L = 0.265 m, e = 9%, / = 90 Hz) 
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FIG. 3. Relaxational dynamics of the wave number k(t) 
in response to a sinusoidal [type(ii)] modulation of the drive 
frequency. The drive amplitude a was co-modulated to keep 
the reduced amplitude e constant. The Faraday pattern re- 
sponds with a sinusoidal variation of its instantaneous wave 
number k(t) with an amplitude Ak m and a temporal phase 
lag 8 (Parameters: L = 0.265 m, s = 9%, / = 90 Hz). 
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FIG. 4. Frequency "dispersion" of the real and imaginary 
parts of Ak* as obtained from modulation experiments of type 
(ii) for the two different container diameters L\ = 265 mm 
and L 3 = 125 mm (/ = 90 Hz, k = 1480 m _1 , e = 9%). 
Solid curves are fits according to the Debye relaxator 

The dependence of the relaxation time r on, respec- 
tively, ho, L, and e can be understood in terms of the 
phase diffusion approach B. Here one takes advan- 
tage of the fact that local disturbances of the elevation 
amplitude die out rapidly, while long-wavelength phase 
perturbation survive on a much longer (diffusive) time 
scale. To streamline the arguments and to work out 
the basic physics we consider a 1-dimensional surface el- 
evation profile in form of stripes as given by Q{x, t) oc 
[ e ife z+¥>(*,i) + c-a ] cos ( w f/2). Here d x <p = Ak describes 
the spatio-temporal variation of the local wave number 
around the underlying base pattern with the wave num- 
ber fco ■ Following the phase diffusion approach the phase 
tp obeys a diffusion equation of the form dtf = D\\d xx ip 
with the diffusion constant (valid close above onset) of 
the form 
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The coefficients t^ 1 = d\/de\k=k c ,e=o and £g = 
— To/2(d 2 \/dk 2 )\k=k a ,e=o are given in terms of the lin- 
ear growth rate A = A(e, fe) at which plane wave per- 
turbations grow out of the plane undeformed interface, 
when the Faraday instability sets in. Here k c is the 
wave number at onset of the Faraday instability. For 
weakly damped (vKq/u) <C 1) capillary waves [kg ~ 
(p/er) 1 / 3 (w/2) 2 / 3 ] on a deep (kh 3> 1) fluid layer (reason- 
able approximations in our experiment) one obtains ap- 
proximately r^ 1 ~ 2vkl and£g = (1/2) [9a / (16v 2 p)]k^ 3 . 
By decomposition of the phase perturbations into a set 
of discrete Fourier modes compatible with the finite con- 
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tainer dimension, ip = X^^Li a n sm i n7T /L), the mode 
n = 1 has the slowest decay time 
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and thus determines the relaxation time of the wave num- 
ber of the experiment. The dashed line in Fig. || is the 
prediction according to Eq. (Q). For a quantitatively 
more reliable check we also evaluated the coefficients To 
and £q numerically from a linear analysis, which takes 
into account the finite filling level as well as gravitational 
contributions to the wave dispersion. The respective re- 
sult is shown by the solid line in Fig. |5|a. Furthermore, 
the predicted quadratic dependence of r on the container 
dimension L is verified by Fig. ^|b. The theoretical pre- 
diction Eq. (||) also implies a slight dependence of r on 
the drive amplitude s. However, checking for this feature 
requires to account for the fact that also the mean wave 
number kg is affected by the drive strength. We deduced 
the empiric dependence ko/k c = 1 — /3s with f3 = 0.264 
from a control run where e was slowly ramped at fixed /. 
Inserting this result into the last term on the right hand 
side of Eq. (0) leads to the following expression 
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Although this relation (see solid line in Fig. ||c) gives 
a reasonable estimate for the e-dependence of r it does 
not correctly reflect the empiric dependence. Apparently 
this is a finite size effect, which is expected to become 
significant at small values of e. Roth et al. (b|] recently 
demonstrated that the effective phase diffusion constant, 
measurable in RBC and TCF relaxation experiments, de- 
pends sensitively on the aspect ratio a = yfeL/^Q, defined 
by the quotient between the container dimension and the 
linear correlation length. Taking TCF as an example, a 
decrease of a from 50 down to 15 (which in our exper- 
iment corresponds to a reduction of e form 9% to 2%) 
implies a decay of the relaxation time by 20 — 30%. This 
is of the same order of magnitude as the value observed 
in our measurements (see Fig. |^c). 
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FIG. 5. The characteristic relaxation time r as a function 
of (a) the mean wave number fco, (b) the container size L, and 
(c) the reduced drive strength e. The dashed line shows the 
prediction for infinite depth capillary waves [Eq. (g)], while 
the solid line is based on a numerical evaluation for the (more 
realistic) case of finite depth gravity-capillary waves. (Param- 
eters: L = 0.265 m e = 9%, / = 90 Hz) 
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